Ultra-cold WIMPs: relics of non-standard pre-BBN cosmologies 
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Weakly interacting massive particles (WIMPs) are one of very few probes of cosmology before Big 
Bang nucleosynthesis (BBN). We point out that in scenarios in which the Universe evolves in a non- 
standard manner during and after WIMP kinetic decoupling, the horizon mass scale at decoupling 
can be smaller and the dark matter WIMPs can be colder than in standard cosmology. This would 
lead to much smaller first objects in hierarchical structure formation. In low reheating temperature 
scenarios the effect may be large enough as to noticeably enhance indirect detection signals in 
GLAST and other detectors, by up to two orders of magnitude. 
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In supersymmetric models, the lightest supersymmet- 
ric particle, usually a neutralino x, is a good cold dark 
matter candidate. This is an example of a more general 
kind of cold dark matter candidates, weakly interacting 
massive particles (WIMPs). The relic density and ve- 
locity distribution of WIMPs before structure formation 
depend on the characteristics of the Universe (expansion 
rate, composition, etc.) before Big Bang Nucleosynthesis 
(BBN), i.e. at temperatures above T ~ 1 MeV. This is an 
epoch from which we have no data. Indeed, if dark mat- 
ter (DM) WIMPs are ever found, they would be the first 
relics from that epoch that could be studied. Signatures 
of a non-standard pre-BBN cosmology that WIMPs may 
provide are few. Here we present one of, to our knowl- 
edge, only three. 

In the standard scenario of WIMP decoupling, one as- 
sumes that the entropy of matter and radiation is con- 
served, that WIMPs are produced thermally, and, of 
most relevance for this work, that the first temperature 
of the radiation dominated epoch before BBN is high 
enough for WIMPs to have reached kinetic and chemi- 
cal equilibrium before they decouple. WIMPs decouple 
first chemically and then kinetically. The chemical decou- 
pling (or freeze-out) temperature T{ a is the temperature 
after which their number practically does not change, 
and in the standard case Tf _ st d — ^x/^O, where m x 
is the WIMP mass. The kinetic decoupling temperature 
Tkd is the temperature after which WIMPs do not ex- 
change momentum efficiently with the cosmic radiation 
fluid. Within the standard cosmology (SC), Tkd-std lies 
between 10 MeV and a few GeV [J. 

There are, however, well motivated cosmological mod- 
els in which the standard assumptions above do not 
hold. These non-standard models include models with 
gravitino Q or moduli [|[ decay, Q-ball decay 0], ther- 
mal inflation the Brans-Dicke- Jordan Q cosmolog- 
ical model, models with anisotropic expansion [jj or 
quintessence domination Q . It has been pointed out that 
in all of these models the neutralino relic density f£v may 
differ from its standard value f2 st d (see e.g. Ref. [9(). 

One clear signature of a non-standard cosmology be- 
fore BBN would be WIMPs that compose part or all 



of the DM but would be overabundant in the SC. Also 
the relic velocity distribution before structure formation 
in the Universe may differ from that in the SC. It has 
already been mentioned in the literature [ToL fTl| that 
WIMP's could be "hotter" than in the SC, even consti- 
tuting warm instead of cold DM. This would leave an 
imprint on the large scale structure spectrum. 

Here we point out a third possible signature of non- 
standard pre-BBN cosmologies: WIMPs may be "colder" 
(i.e. they may have smaller typical velocities and thus 
smaller free-streaming length) and the mass contained 
within the horizon at kinetic decoupling may be smaller 
than in the SC. This would lead to a smaller mass for the 
smallest WIMP structures, those formed first. Some of 
the smallest WIMPs clumps would survive to the present. 
Smaller and more abundant DM clumps would be present 
within our galaxy, an observable consequence of which 
would be a stronger annihilation signal from our galac- 
tic halo detected in indirect DM searches by GLAST, 
PAMELA and other experiments [H, [H [3 111 El, E3, 
UH. The signal in direct DM searches might also be af- 
fected in significant ways [l!|. We show that in low re- 
heating temperature cosmological models, both the free- 
streaming mass scale and the mass within the horizon 
at kinetic decoupling may be smaller than in the SC by 
many orders of magnitude 

We present estimates of the kinetic decoupling temper- 
ature, characteristic relic WIMP velocity, free-streaming 
and kinetic decoupling horizon mass scales using a 
generic WIMP elastic scattering cross section written as 
in Ref. [14[ and order of magnitude calculations. Af- 
ter chemical decoupling, T < Tf D , the total number of 
WIMPs remains constant and WIMPs are kept in local 
thermal equilibrium by elastic scattering with relativistic 
particles in the plasma. The WIMPs are non-relativistic 
at these temperatures, thus the average momentum ex- 
changed per collision is small, of order T, and the rate of 
momentum exchange T is suppressed by a factor ~ T/m x 
with respect to the rate of elastic scattering 
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Here a e \ is the total cross section for elastic scattering of 
WIMPs and relativistic Standard Model fermions, n rad ~ 
T 3 is the number density of relativistic particles, which 
are assumed to be in local thermal equilibrium, and v ~ 1 
is the WIMP-fermion relative velocity. The thermal aver- 
age of Cei can be written as (<7 e i) = <jf(T/m x ) 1+l 1 where 
erg 1 ~ {Gym^)' 1 rr? Irn\ ~ 10~ 10 m 2 GeV~ 4 sets the mag- 
nitude of the cross section, and I parametrizes its temper- 
ature dependence. Finally, m\y and mz are the masses 
of the standard gauge bosons and Gp is Fermi's coupling 
constant. In the Standard Model, elastic scattering be- 
tween a light fermion and a heavy fermion is mediated by 
Z exchange and 1 = 0. In supersymmetric extensions of 
the Standard Model, where the lightest neutralino is the 
WIMP candidate, sfermion exchange occurs if the neu- 
tralino is a gaugino, Z exchange is suppressed, and 1 = 1. 
Inserting erg 1 into Eq. [TJ the rate of momentum exchange 
for non-relativistic WIMPs is 

10-iogB / T y 

" GeV 4 \m x ) ■ (2) 

In the following, we focus on the case of neutralinos with 
I = 1 (analogous results can be easily derived for 1 = 0). 

Kinetic decoupling occurs when the rate of momentum 
exchange becomes smaller than the expansion rate of the 
Universe H. In the SC, decoupling happens while the 
universe is radiation dominated so the Hubble parameter 
is H ~ T 2 /M P , where M P ~ 10 19 GeV is the Planck 
scale. From r ~ H, we get 

T kd _ std .20Me V ( I ^ v ) . (3) 

More accurate calculations give a range of 10 MeV to a 
few GeV for T kd _ std Q. 

We concentrate now on a class of non-standard cosmo- 
logical models with a late episode of inflation or entropy 
production 0, 0, H, 0| in which a scalar field (ft dominates 
the energy density of the Universe and subsequently de- 
cays (while oscillating around a minimum of its potential) 
eventually reheating the Universe to a low reheating tem- 
perature Trh- This does not s poil primordial nucleosyn- 
thesis provided Trh ^4 MeV [20]. The interesting case 
for us is when Trh is smaller than the standard chemi- 
cal decoupling temperature Tf _ s td, so kinetic decoupling 
happens during the c6-oscillations dominated phase. 

Late-decaying scalar field models are well motivated 
in particle theories. For example, moduli fields which ac- 
quire a mass at the supersymmetry breaking scale 10 
to 100 TeV and have gravitational strength interactions, 
thus their decay rate is r docay ~ to^/M 2 , are pervasive 
in supersymmetric models. These fields naturally tend to 
dominate the energy density of the Universe at late times 
and produce reheating temperatures in the MeV range 
(the "moduli problem" [2l[ is the tendency of these de- 
cays to happen even after BBN, which must be avoided). 
In fact, approximating the decay as instantaneous (usu- 
ally a very good approximation) at the moment of decay 



the energy stored in the field goes into radiation at a 
temperature T RH . Thus r docay ~ -ff(T RH ) ^ T RH /M P 
emplies that 

T ™^ M 'Hv£hsT 2 - < 4 > 

During the epoch in which the Universe is dominated 
by the oscillating (f> field, the Hubble parameter H$ is 
proportional to T 4 (22j. Since at the moment of <f> de- 
cay, when T = T RH , H^Trb.) — ?r H /Mp, we can fix 
the proportionality constant so that ~ T 4 / (T^ H Mp). 
Requiring that T ~ at the new kinetic decoupling 
temperature Tkd', we obtain 

Combining the two equations r(Tk d - s td) — T^ d _ std /Mp 
and r(T kd ,) ~ T* d ,/{Tl H M P ), we have 

T 2 

rri kd — Std / n\ 

T kd ' ■ (6) 

Thus, if the reheating temperature is smaller than 
the standard kinetic decoupling temperature (T R h < 
Ikd-std), WIMPs decouple earlier than in the SC, i.e. 
Tkd' > Tk d -std, and do so during the c6-oscillations dom- 
inated epoch. 

So far we have assumed that the WIMPs are non- 
relativistic at decoupling. Thus, the relations above hold 
for Tk d ' < m x /3. If Eq. [5] leads to Tk d ' > m x /3, the 
WIMPs would be relativistic at decoupling, and the equa- 
tions need to be modified. Since the momentum trans- 
fer to the radiation background is very efficient in col- 
lisions of relativistic WIMPs, in general the kinetic de- 
coupling would happen at the moment WIMPs become 
non-relativistic and not earlier, i.e. Tk d ' ~ m x/^ (unless 
the scattering cross section is so small that WIMPs are 
never in kinetic equilibrium [11|). 

Late-decaying scalar field models exemplify many com- 
binations of reheating and decoupling temperatures. In 
these models, the dominant WIMP production mecha- 
nism can be thermal (due to interactions with the ra- 
diation background) or non-thermal (due to the decay 
of the 4> field into WIMPs), with or without chemical 
equilibrium (see for example Ref. @). In these models, 
neutralinos in almost all supersymmetric models could 
have the relic density necessary to be the DM Q through 
a combination of thermal and non-thermal production 
mechanisms 0, H, H H H M, H, M, HI- For 
thermal production without chemical equilibrium, most 
WIMP production happens at T* ~ m x /4 [25j]. Thus 
the WIMP number per comoving volume is fixed then. 
For Tk d ' < Tt, the equations we derived above hold. For 
thermal production with chemical equilibrium the neu- 
tralino freezes out while the Universe is dominated by 
the 4> field at a new freeze-out temperature Tf Q > higher 
than the usual m x /20 [22|, [26j]. The freeze-out density is 
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FIG. 1: Ratios of the late decaying scalar field (with ul- 
tra cold (UC) WIMPs) and standard (std) scenario free- 
streaming scales (Mf£ c /Mj£ d , solid lines) and damping scales 
(M^ c /Mf d , dashed lines), as functions of the standard ki- 
netic decoupling temperature Tkd-std- Each set of two lines is 
labeled by the corresponding value of the reheating tempera- 
ture Trh- 



larger than usual, but it is diluted by entropy produc- 
tion from cj) decays, fi x ~ T^ H Tf _ s td(T f( y )~ 4 f2 st d- The 
numerical results in Ref. 9] indicate a dependence of fl x 
closer to Tpj H , and that the freeze-out temperature de- 
pends of T RH as T io >20/m x ~ {T RH 20/m x )^ 4 . In this 
case, the equations we derived hold for Tkd' < 2fo/. Non- 
thermal production without chemical equilibrium hap- 
pens when production of WIMPs in the decay of the <j> 
field is not compensated by annihilation. WIMPs are 
produced with an energy which is a fraction / of the <f>- 
field mass m^, E x ~ f m^. Thus if m x < 3/to0, WIMPs 
are produced relativistic. In this case the elastic scatter- 
ing cross section is (va c \) — aoTE x /m x ~ aoTfm^/rn^. 
Taking the characteristic value of cro as above, neutrali- 
nos are in kinetic equilibrium while they are relativistic 
for m > 10 keV//(T RH /10 MeV) 2 , that is for all the 
physically acceptable values of m§. Thus kinetic decou- 
pling occurs while neutralinos are non-relativistic as as- 
sumed above. In the extreme case in which neutralinos 
are never in kinetic equilibrium [Tl| . neutralinos can ac- 
tually be warm DM 10], since they are produced late in 
the history of the Universe and with a large initial energy 
that redshifts until the moment of structure formation. 

We now estimate the WIMP characteristic speed in 
low Trh models. At the moment of kinetic decoupling, 
WIMPs are in thermal equilibrium with the radiation, 
and their characteristic speed is v'(T kd ') — y/T^ /m x . 
After decoupling, the speed decreases as a -1 (if the 



WIMPs are non-relativistic). During the 0-oscillations 
dominated era, the scale factor of the Universe is related 
to the radiation temperature as a ~ T"®^. When the 
Universe becomes radiation dominated, i.e. at Trh, the 
characteristic speed of WIMPs is therefore 



w'(Trh) ^ \ T kd ,/m x (Trh /Tkd') 



8/3 



(7) 



In comparison, in the standard radiation dominated case, 
a ~ T _1 and at the same temperature Trh WIMPs have 
a characteristic speed 

v s td(TRn) ~ \jT kd - std /m x (TR H /T k d- S td) • (8) 

Because speeds redshift in the same way in both models 
at temperatures smaller than Trh, after reheating the 
speeds v' and v s td remain in the same ratio 



v' /v std ~ (TR H /T kd -std) 



10/3 



(9) 



This relation applies to the case Ted' < v ^xl^ > f° r wu ich 
Eq.[S]holds. Thus the characteristic relic WIMP speed in 
low Trh cosmological models can be much smaller than 
in the SC. In other words, WIMPs can be much colder, 
i.e. "ultra-cold" , as we call them. 

The free-streaming length Af s of ultra-cold WIMPs is 
consequently smaller than that of standard WIMPs. Af s 
is the characteristic distance covered by WIMPs from the 
time of kinetic decoupling tkd to the present (while they 
propagate as free particles) 



[ t0 dt 
A fs = ca I v- 



tkd 




a 3 H{a) ' 



(10) 



Here, Tkd and GSkd are the temperature and scale factor 
at the moment of kinetic decoupling, and ao is the value 
of the scale factor at present. During the ^-oscillations 
dominated epoch, H(a) oc a~ 3 / 2 with a ~ t 1 ! 3 ~ j"-8/3 
and H ~ T A / (T^Mp). During the radiation dominated 
epoch, H(a) oc a" 2 , with a ~ i 1/2 ~ T" 1 and H ~ 
T 2 /Mp, and the integral in the definition of Af s is oc 
In a. During the matter-dominated epoch, H (a) oc a~ 3 / 2 , 
with a - i 2 / 3 ~ T -1 and H ~ T 3 / 2 /Af P , and the free- 
streaming length saturates. 

When the kinetic decoupling occurs during the cf>- 
oscillations dominated epoch and WIMPs are ultra-cold 
(UC), we obtain 
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1.96^. (11) 



The first term within the curly brackets arises from the 
(^-dominated epoch, the logarithmic term from the ra- 
diation dominated epoch, and the term 1.96 = 2[1 



(1 



J 4 / 2 ] from the matter dominated epoch. The 



4 



subindex eq refers to matter-radiation equality. When 
Eq.HJholds (for T kd , < m x /3), Eq. [TT] becomes 
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cM T 



RH 



t kd-std 
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Tkd-std 
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Id [ 5^ + 1.96 i .(12) 
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This is to be compared with the free-streaming length in 
the SC, 
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cM r 



7o^/m x Tkd-std 



In 



kd-std 



T, 
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1.96 



(13) 



As traditional, we introduce the mass Mf s contained 
within a sphere of radius Af s /2, and compare the stan- 
dard and non-standard scenarios through the ratio 



M std 



\UC 
A fs 

A fs 



(14) 



As shown in Fig. 1 (solid lines), this ratio is always 
smaller than 1 if the reheating temperature is smaller 
than the standard kinetic decoupling temperature, and it 
can be many orders of magnitude smaller. For Trh = 5 
MeV and T kd _ std between 10 MeV and a few GeV Q, 
the free-streaming scale Mf s can decrease by a factor be- 
tween 0.1 and 10 -13 (the suppression is less important 
for larger values of Trh)- 

Friction between WIMPs and relativistic leptons dur- 
ing kinetic decoupling (Silk damping) leads to a small- 
scale cutoff in structure formation at the scale of the 
horizon at kinetic decoupling (l5l . fl6| . The mass con- 
tained within the horizon at decoup ling in the SC is 
Mf d ~ 1O- 4 M (1O MeV/T kd _ std ) 3 [HI [H|. It varies 
from 1O~ 4 M for T kd _ std ~ 10 MeV to 1O" 12 M for 
T k d- S td — 5 GeV . 

During the ^-oscillations dominated phase the Uni- 
verse expands, and thus the density contrast of DM inho- 
mogeneities grow, in the same way as in a matter domi- 
nated phase. A detailed study (which is beyond the scope 
of this paper) of the kinetic decoupling during this phase 
should be done to find the cut-off mass scale of the small- 
est dark matter structures, which will also depend on the 
particular particle physics model considered. However, it 
is reasonable to assume that also in this case the cut-off 
will be given by the comoving free-streaming mass scale 
and/or the kinetic decoupling horizon mass scale. 

During the (^-oscillations dominated phase the mat- 
ter density scales as T~ s and the time as T~ 4 . Thus, 
the mass contained in the horizon at decoupling M^ c is 
smaller than in the SC by the factor 



M std 



?kd' 



kd-std 



(15) 



Using Eq. [5J this ratio becomes (?RH/?kd-std) 5 - As 
seen in Fig. 1 (dashed lines), the suppression factor 



of M} 10 is now from 10" 27 Af Q for T kd _ std 
10- 5 M Q for T kd _ 



(My c /Mj td ) can be substantial, ranging from 0.1 for 
T kd _ std = 10 MeV to 10~ 15 for T kd _ std = 5 GeV, when 
assuming Trh = 5 MeV (the suppression is less impor- 
tant for larger values of Trh)- This means that the range 

5 GeV to 

kd _ std ~ 10 MeV. 
From the ratios (M^ c /Mf d ) and (M^ c /Mg d ) shown 
in Fig. 1, we see that the damping mass scale Md is less 
suppressed than the free-streaming mass scale Mf s , ex- 
cept possibly for standard kinetic decoupling tempera- 
tures above 1 GeV. This is important because the scale 
of the smallest WIMP haloes will be the larger of Mf s and 
Md, since a halo mass must be larger than both. In the 
standard cosmological scenario, the damping scale MJ td 
is usually larger than the free-streaming scale (by a fac- 
tor (m x /T kd _ std ) 3 / 2 , at least for supersymmetric gaugino 
models with T kd _ std in the 10-100 MeV range [lEEl)- 
Fig. 1 shows that in low Trh models the damping scale 
remains larger than the free-streaming scale for most typ- 
ical values of the standard kinetic decoupling tempera- 
ture T kd _ std <1 GeV. Thus, we take the damping scale 
Md to be the characteristic mass of the smallest WIMP 
halos. 

If the smallest halos survive until today, they may 
be present in the dark halo of our galaxy and may en- 
hance the expected WIMP annihilation signals over the 
smooth halo expectation by a boost factor B. The B 
factor increases slowly with decreasing Md- For a halo 
of mass M and smallest subhalo mass Md, Ref. [3lT ] 
finds B ~ 0.1[(Af/M d ) 013 - 1]. In the SC, one ex- 
pects B of the order of 10 for the Milky Way for which 
M ~ 1O 12 M . For example, from the equation just men- 
tioned we get B ~ 20 to 130 for the standard range of 
M d , from 10- 6 Af o to M d ~ 1O~ 12 M . With ultra-cold 
WIMPs, Md could be much smaller and thus the boost 
factor could be much larger: it could reach B ~ 10 4 for 
Md — 1O~ 27 M . Such large boost factors would not only 
make a halo WIMP annihilation signal easier to detect, 
but would also be a signature of a non-standard pre-BBN 
cosmology. 

We finally remark that ultra-cold WIMPs may arise in 
all models in which the expansion rate of the Universe 
in the pre-BBN era is larger than assumed in the SC, 
although the magnitude of the effect would in general 
be smaller than in the low reheating temperature models 
presented above. For example, let us consider "kination" 
models Q. These are models in which the kinetic en- 
ergy of a scalar field, = <fi 2 /2 ~ a~ 6 , dominates the 
energy density of the Universe at T > T k ; n before BBN, 
while the entropy is dominated by the radiation, thus 
a - T _1 . Roughly, p^/p^a ^ ^(T/MeV) 2 , where ^ 
is the fraction of the energy density of the Universe at 
T ~ 1 MeV due to the kinetic energy of the <f> field. 
At higher temperatures the fraction of <f> kinetic energy 
grows very fast with the temperature and it is dominant 

— 1/2 

for T > T k j n ~ 7]^ ' MeV. If the kinetic decoupling 
of WIMPs happens during the kination period, i.e. if 
7k d - s t d > Tkin, assuming H ~ ^fp^/Mp, the kinetic de- 
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coupling happens approximately at 

r^SOMeV,^^) 1 ". (16) 
In this case the free-streaming length is 

fs V m * «oT 2 I T k k - 

ln(^)+1.96|,(17) 

which for T^} > Tkd-std > 2\j n is smaller than the 
standard free-streaming length. During the kination pe- 
riod, the scale factor is a ~ and the time evolves as 
t ~ T~ 3 thus the mass contained within the horizon at 
kinetic decoupling Mf n , again for Tgp > T kd _ std > T kln , 
is smaller than the standard mass scale M^ td by the ratio 

1U d ^ -'kin^kd-std (18) 

Thus in kination models the free-streaming and damping 
mass scales can be smaller than in the SC, although not 
by as much as in the late decaying scalar field (or low 
reheating temperature) models presented above. 

We have here pointed out that a too large boost factor 
in the annihilation signal of a particular WIMP would be 



a signature of a non standard cosmological evolution of 
the Universe just before BBN, during the kinetic decou- 
pling of WIMPs. If dark matter WIMPs are ever found, 
they would be the first relics from the pre-BBN epoch 
that could be studied. Signatures of a non-standard pre- 
BBN cosmology that WIMPs may provide are few and 
here we presented one of them: WIMPs may be ultra 
cold so the mass of the smallest WIMP structures, those 
formed first, may be smaller than in the standard cosmol- 
ogy. Some of the smallest WIMPs clumps would survive 
to the present. Smaller and more abundant DM clumps 
would be present within our galaxy, an observable conse- 
quence of which would be a stronger annihilation signal 
from our galactic halo detected in indirect DM searches 
by GLAST, PAMELA and other experiments. Boost fac- 
tors as large as 10 4 for usual WIMP candidates are pos- 
sible in the low reheating temperature scenarios consid- 
ered here. In last instance, verifying this signature would 
require to study in accelerators, the LHC or ILC, the 
properties of the particular WIMP that would allow us 
to estimate its scattering cross section, to find the same 
WIMP in indirect detection searches and to understand 
the formation and survival of the earliest dark matter 
clumps within the halo of our galaxy. 
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